Unconventional magnetism in small gold organic molecules 
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We present a theoretical study of the magnetic properties of dicyclopentadienyl metallocene and phthalo- 
cyanine molecules, that contain the transition metal atoms M = Fe, Co, Ni, Cu, Zn, Ir, Pt and Au. Our most 
important prediction is that gold and copper molecules are magnetic. We find that the magnetism of these 
molecules is fairly unconventional: the gold atom itself is weakly magnetic or even non-magnetic. Its role is 
rather to induce magnetism in the surrounding carbon and nitrogen atoms, producing a sort of spin density wave. 

PACS numbers: 75.75.+a, 33.15.Kr, 81.07.Nb 



The efforts devoted in the past few years to fabricate and 
characterize new nano- scale objects have helped to uncover a 
wealth of fascinating geometric, mechanical, electronic, mag- 
netic, optical or dissipative properties, that are brought about 
in many cases by the laws of quantum physics. More specif- 
ically, the multidisciplinary field of Nanomagnetism aims 
at the fabrication of devices with tailored magnetic proper- 
ties. Indeen, there exists an intense activity in the micro- 
magnetic community to extend its reach to the nanoscale, by 
reducing the size of their room-temperature magnetic dots 
well below the 100 nm mark. Gold nanoparticles of diam- 
eters in the 2 nanometer range have been found recently to 
show room-temperature magnetism, when capped with or- 
ganic moleculestll|2|,[3|]. Interestingly, the experiments have 
found that the magnetism of these nanoparticles is localized at 
their surface, where Au atoms are in chemical contact with the 
capping molecules. The different experiments produce differ- 
ent estimates for the spin moments Ms of the surface gold 
atoms, which vary from 0.002 to 0.3 /isis)], but they agree on 
the very high value of the magnetic anisotropy energy (MAE), 
which is of about 0.4 eV/atom. Huge orbital moments Ml at 
the surface gold atoms have been proposed to exist and orig- 
inate that large MAE|0], but the actual measured moments 
seem to be modest, of order Ml /Ms ~ 0.15|[ll|3|]. 

Molecular magnets are in some ways superior to nanopar- 
ticles, since they have a well defined number of atoms with 
precise chemical identities, and do not suffer from particle 
number dispersion. Further, molecules are not so prone to 
conformational changes, since these require breaking a siz- 
able number of covalent bonds, as opposed to the metallic 
bonds of atomic clusters. Finally, the chemical activity of 
molecules can be engineered by oxidation or reduction, or by 
the attachment of end groups. Molecular magnets contain- 
ing single Rare-Earth ions have been studied in the pastJUB], 
but late 5d elements like iridium, platinum or gold are pos- 
sibly better candidates for room-temperature magnetism since 
these are heavier elements, which means that they have higher 
Spin-Orbit coupling constants. Our previous calculations for 
platinum dimers found MAEs of the order of 0.1 eV/atomj7|], 
which are consistent with the experimental values in gold 
nanoparticles referred above. 

We present here a detailed theoretical study of the mag- 



netism of two popular organic molecules, metallocenes 
and phthalocyanines. Metallocenes, denoted MCp2, are 
organometallic molecules that contain two cyclopentadyenil 
rings Cp, which sandwich a transition metal atom M. Metal- 
lic phthalocyanines, denoted MPc, are macrocyclic molecules 
that have an alternating nitrogen atom-carbon atom ring struc- 
ture, and contain a metal ion M in its center, which bonds with 
the four isoindole nitrogen atoms. We have used the third row 
elements Fe, Co, Ni, Cu and Zn, and the fifth row elements 
Ir, Pt and Au as the metallic atom M. The axial symmetry of 
these molecules produces a crystal field at the metal ion posi- 
tion which leads to expect that Fe and Zn molecules have a net 
zero spin, that Co, Ir, Cu and Au molecules have spin 1/2, and 
that Co and Pt molecules have spin 1 . Our simulations con- 
firm these guesses except for Pt molecules, that have a spin- 
less ground state. Interestingly, we find that gold molecules 
have a finite spin 1/2, but that the distribution of magnetic 
moments across gold molecules is not conventional. The gold 
atoms produce a sort of spin density wave such that the sur- 
rounding carbon or nitrogen atoms are also magnetized, show- 
ing a spin-density-wave pattern, with almost collinear orienta- 
tions. The gold atom itself contributes about 15 % of the total 
spin moment of the molecule in orocene, while it is unmagne- 
tized in the phthalocyanine. The orbital moments are between 
three and thirty times smaller that the spin moments. We have 
checked that all the spin 1/2 molecules show zero magnetic 
anisotropy as they should lll3ll . The spin-1 Ni molecules also 
show very small anisotropies, which was expected because of 
its small Spin Orbit coupling constant. 

The structural and electronic properties of metallocene and 
possibly also of phthalocyanine molecules in the gas and crys- 
talline phases are well captured by density functional the- 
ory (DFT)lIi. We have therefore carried the simulations 
with the molecular dynamics SIESTA suite(9|], which uses 
norm-conserving pseudopotentials and a basis set of localized 
atomic-like wave functions. We have used a triple-zeta doubly 
polarized basis set for transition metal ions (Fe, Co, Ni, Cu, 
Zn, Ir, Pt and Au), and a double-zeta polarized basis set for the 
carbon, nitrogen and hydrogen atoms. We have checked ex- 
plicitly that our simulation boxes for metallocenes, having 20 
A of lateral size, where large enough to avoid spurious bound- 
ary effects of the electric field. For the larger phthalocyanine 
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FIG. 1 : (Color online) (a) eclipsed, and (b) staggered, geometries of 
metallocene molecules, showing rf coordination; (c) phthalocyanine 
molecule. 



molecules, the boxes were cuboids of size 28x28x20 A 3 . We 
have used a very fine grid for the real space integrals to ensure 
that egg-box effects were absent. We found it essential to set a 
rather strict tolerance of 5 x 10 -4 eV/A in the force relaxation 
procedure of the atoms, in order to scape from local minima of 
the energy landscape. We took between 5 and 10 different ini- 
tial seeds for the geometry /spin arrangements in our sampling 
of the Hilbert space of each quantum system to ensure that 
we would reach the ground state configuration. We have con- 
firmed that the Local Density and the Generalized Gradient 
approximations (LDA and GGA, respectively) 03] shed 
the same results for the ground state and first excited isomers 
of all the molecules, with only small quantitative differences. 
We have also double-checked our results for gold molecules 
with the alternative code QUANTUM ESPRESSOlHl, which 
uses ultra-soft pseudopotentials and a plane- wave basis set. 

Figs. 1(a) and 1(b) show the eclipsed and staggered ge- 
ometrical structures of metallocene molecules, which have 
the standard rf coordination, whereby the central metal ion 
bonds with all five carbon atoms in each ring. There exist 
other possible geometrical arrangements, where the rings are 
displaced laterally, or even tilted, so that only one, two or 
three carbons in each ring bond to the metal ion; the coordina- 
tion is called 77 1 ' 2 ' 3 , accordingly. All of the metallocenes that 
we have simulated have been synthesized, except for AuCp2, 
and possibly CuCp2 . Their geometrical structures for the gas 
phase have been determined experimentally for Fe, Co, Ni 
and Zn|[3, 15]. Detailed quantum chemistry simulations of 
the geometry and total spin of isolated Fe, Co, Ni, Cu and Zn 
metallocenes have been performedll8[ 1 1 6l 1 1 7l 1 1 811 . which agree 
with the experimental data. Our results both for the geometry 
and the spin ground state, which are summarized in Table I, 
agree in detail with the available experimental and theoretical 
data. 

We select the cobalt, copper and gold metallocenes for 
a closer inspection of the magnetic phenomena of these 
molecules. The rationale behind our choice is the following. 
First, the three molecules have the same spin (S=l/2) because 
of the axial crystal field of the molecule, so that their magnetic 
properties can be meaningfully compared. Second, cobalt is a 
strong ferromagnet, while copper and gold have a noble metal 
electronic shell structure. Third, gold has more delocalized s- 
wave functions, so it is expected to show a weaker magnetism. 



Molecule 


Exp. geometry 


Our geometry 


(an) 




FeCp2 


Eclipsed rj b 


Agree 


0.0 


0.0 


CoCp 2 


T- i; 1 1,2 

Eclipsed 77 ' 


Agree 


0.99 


0.75 


NiCp2 


Eclipsed 77 5 


Agree 


1.99 


1.22 


CuCp2 


Not synthezised 


Staggered 77 2 — 77 5 


f\ no 


U.ZZ 


ZnCp2 


Staggered 77 1 — 77 5 


Agree 


0.0 


0.0 


IrCp 2 


Not available 


Eclipsed 


0.95 


0.50 


PtCp 2 


Not available 


Staggered 77^ — rf 


0.0 


0.0 


AuCp 2 


Not synthezised 


Staggered 77 1 — rf 


0.99 


0.13 



TABLE I: Structural and magnetic properties of the metallocene 
molecules simulated in this article. The total spin moment of the 
molecules is denoted by Mt and is oriented along the z-axis. The 
spin moment of the transition metal atom Mm is estimated from a 
Mulliken population analysis. 



Fourth, we wish to understand better the magnetic properties 
of small gold molecules in view of the recent experiments 
mentioned abovelH 0, 0]. We plot in Fig. 2(a) a schematic 
view of the geometrical structure and atomic spin moments of 
CoCp2. We have found that the most stable isomer has basi- 
cally D$h symmetry, with small Jahn-Teller distortions, which 
favor ?7 2 coordination. More important for the present discus- 
sion is the crude size of the atomic moments in the molecule, 
which the figure shows to be localized in the cobalt atom. Ta- 
ble I indeed shows that Co carries 75% of the total moment 
of the molecule (Mc = 0.75 lib)- Notice that the rf — rf 
coordination is reflected in the small moment of two of the 
carbon atoms in each carbon ring. One of the two atom rings 
in the D$d geometry CuCp2 has been displaced and shows rf\ 
the other ring shows a much more slight displacement, which 
leads to rf coordination, as we show in Fig. 2(b). Further- 
more, the two rings are not coplanar. The magnetic configura- 
tion of this molecule is rather interesting. The transition metal 
atom has a much lower magnetic moment, which amounts to 
only 20 % of the total moment of the molecule, see Table I. 
Fig 2(b) shows that the magnetization of the carbon atoms 
follows a non-monotonic dependence with the distance to the 
copper atom, which is reminiscent of a sort of spin-density 
wave. Indeed, two of the carbon atoms in the displaced ring 
have spin moments of 0.25 lib, while the other three are very 
weakly magnetized (M ~ 0.02/i^). A similar situation hap- 
pens in the other ring, where two carbon atoms have sizable 
magnetizations of about 0.13 lib, while the moments of the 
other three are tiny. This unconventional behavior is more 
apparent in AuCp2, that we show in Fig. 2(c). The two cy- 
clopentadyenil rings are displaced in opposite directions by 
the same distance, which leads to rf — rf coordination. The 
two rings are coplanar in this case, which makes a much more 
symmetric molecule than CuCp2 . The figure clearly indicates 
that two carbon atoms in each ring have a much larger mag- 
netization (0.21 lib) than the one shown by the central gold 
atom (0.13 lib)- The magnetization profile shows a more pro- 
nounced non-monotonic behavior than CuCp2i the two car- 
bon atoms in each ring which are directly bonded to the gold 
atom are unmagnetized, while the carbon atom furthest apart 
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FIG. 2: (Color online) Geometry and atomic spin moments of the 
most stable isomers of (a) CoCp2; (b) CuCp2; and (c) AuCp2. No- 
tice that the energy is unaffected by a global spin rotation, since the 
magnetic anisotropy energy is essentially zero in the three molecules. 
We have chosen to display here a case were the spins are oriented 
perpendicular to the cyclopentadienyl rings. 



shows antiparallel magnetization. This suggest even more 
strongly than in the CuCp2 case that we are facing with a spin- 
density wave profile. 

We note now that we have used a Mulliken population anal- 
ysis to draw the above conclusions on the magnetism of metal- 
locene molecules. To dispel fears on the possible inaccuracies 
of this analysis, we have estimated also the atomic charge and 
spin moments by integrating the charge and spin densities in- 
side spheres centered in each atom, for a number of different 
integration radii, that we show in Fig. 3(a). We find that the 
spin density curves show plateaus for radii slightly larger than 
the inter- atomic distances. The values of the spin moments 
at these plateaus coincide almost quantitatively with the esti- 
mates provided by the Mulliken analysis. The plateaus are not 
present for the atomic charges, apart from a small shoulder in 
the case of gold. To understand better the spatial distribution 
of the integrands, we plot in Figs. 3(b) and (c) the density 
iso-contours. The charge iso-contours show the charge con- 
jugation appropriate to the covalent 7r chemical bond among 
the p z orbitals in the carbon rings. This is the reason why the 
integrated charges of the carbon atoms do not show plateau 
behavior. In contrast the charge about the central gold atom 
has a spherical shape, which points towards a more ionic char- 
acter of the gold-carbon bond. This is also reflected in the 
shoulder referred above. The spin iso-contours have a lobu- 
lar structure. The four lobes that appear for Au correspond 
to the d zy orbital, which is the most directly involved in the 
gold-carbon chemical bond responsible for the rj 2 coordina- 
tion. The lobular structure of the carbon atoms corresponds 
to the p z orbitals; the lobes point slightly towards the ad- 
jacent carbon atoms, which indicates a larger derealization 
due to the conjugation present in the ring. We note that the 
atoms that showed antiparallel spin moments in each ring in 
Fig. 2(c), show here a negative spin density. The localiza- 
tion of the spin density in the orbitals making up the chemical 
bonds in the molecule, together with the non-uniform varia- 
tion of the atomic moments strongly suggest again that the 
spin landscape in this molecule corresponds to a static spin 
density wave originated by the presence of the gold atom. 

We finally note that the magnetic moments are essentially 




FIG. 3: Atomic spin moments and charge of AuCp2. (a) Charge and 
spin moments estimated by integrating the charge and spin densities 
on spherical contours. The abscissa in the graph corresponds to an 
imaginary line that begins at the Au atom, and goes to one of the 
cyclopentadyenil rings, where it passes first through the closest car- 
bon atom(C(l)), then through the following one (C(2)) and ends up 
on the farthest (C(3)). The position of each of these carbon atoms 
is marked by a vertical line in the graph. The colored rectangular 
bars indicate the Mulliken population results for the charge and spin 
moments for each atom. The different data points correspond to in- 
creasing integration radii for the spherical integration contours, (b) 
and (c) Charge and spin iso-contours of the molecule. 




FIG. 4: (Color online) Geometry and atomic spin moments of the 
most stable isomers of (a) CoPc; (b) CuPc; and (c) AuPc. The plane 
of the CoPc and CuPc molecules is tilted in the figure on purpose. 
Notice that the spins in AuPc show small non-collinearities. 



collinear. We have also computed the expectation value of the 
atomic orbital moments, which we find to be 0.02-0.03 \ib 
at the gold atom, and below 0.001 /i^ at most for the sur- 
rounding carbon atoms (as otherwise expected since the spin 
orbit coupling constant is only large for Au). We therefore 
find that the ratio M L /M S is of the order 0.2-0.3 at most. 
This conclusion agrees with the experimental results for gold 
nanoparticles [Q U . 

The above trends on the magnetic behavior of metallocenes 
molecules can be generalized to other organic molecules con- 
taining transition metal ions, like phthalocyanines. We have 
found that the geometry of these molecules does not change 
during the force relaxation cycle, and only the positions of 
the central atoms are slightly modified. Cobalt, copper and 
gold phthalocyanines (CoPc, CuPc and AuPc, respectively) 
all have spin 1/2 and zero magnetic anisotropy, as we desired. 
We have therefore chosen to plot in Fig. 4 our results for 
molecules laying in the XY plane, with spins oriented along 
the z-axis for CoPc and CuPc; for AuPc in contrast, we found 
it more illustrative to plot a case where the spins were oriented 
in the plane of the molecule, and along the x-axis for AuPc. 
We find that the atomic spin moments are very strongly lo- 
calized at the cobalt atom for CoPc, but that they are spread 
towards the neighboring nitrogen atoms in CuPc. Gold ph- 
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FIG. 5: (a) A top view of the spin density wave in one of the two rings of orocene. Red (blue) colour indicates positive (negative) spin densities. 
The colour intensity indicates how positive or negative is the spin density, (b) Top view of the charge iso-contours in a gold phthalocyanine. 
(c) and (d) Lateral and top views of the spin density wave in AuPc. The colour code is the same as in AuCp2. 



thalocyanine is specially interesting again. The spin moments 
are spread at the central core and all along the benzene rings 
in one of the arms of the molecule, while the moments of the 
central gold atom and of the benzene rings at the other arm are 
negligible. Hence, the spin configuration in the ground state 
breaks the C^ v symmetry of the molecule down to C2. We 
stress that the moment at the gold atom is negligible, and that 
the role of this atom seems to produce a spin-density wave, 
magnetizing its surrounding atoms. As in gold metallocene, 
we find that the orientation of the magnetic moments in AuPc 
is collinear to a first approximation. Interestingly, we have 
found slight non-collinearities when we orient the spins in the 
plane of the molecule, as is apparent in Fig. 4(c). In contrast, 
the non-collinearities dissappear when the orientation of the 
atomic spins is perpendicular to that plane. 

We finally wish to pay a closer look at the spin-density pro- 
file of AuCp2 and AuPc. We plot to this end the spin den- 
sity in one of the cyclopentadyenil rings of orocene, in Fig. 
5(a). We find that this density is largely spread throughout the 
ring, as opposed to being strongly localized about the carbon 
atoms, and shows wide regions of opposite charge. We plot 
in Fig. 5(b) the charge density of AuPc, which shows the co- 
valency of the benzene rings, and the more ionic character of 
the bonds linking nitrogen and gold atoms. Fig. 5(c) shows a 
lateral view of the spin density in AuPc, which indicates that 
the magnetization is localized in the p z nitrogen and carbon 
orbitals. A top view of the molecule in Fig. 5(d) indicates that 
the spin density is also delocalized in this molecule, showing 
small regions of negative spin surrounded by much larger re- 
gions where the spin density is positive. 

As a summary, we find that small gold organic molecules 
show unconventional magnetic features. We find that the role 
of the gold atom is to magnetize the surrounding atoms, while 
displaying themselves small or even negligible magnetic mo- 
ments. The magnetization pattern shows a spin-density wave 
behavior, whereby the spin profile is fairly delocalized across 



the molecule, and presents regions with positive and with neg- 
ative spin density. We find that the orientation of the spins 
might be parallel or anti-parallel, but is collinear to a first ap- 
proximation. We also find that the orbital moments are much 
smaller that the spin moments. 

We are indebted to I. K. Schuller and F. J. Garcia Alonso. 
We acknowledge financial support from the Spanish Research 
Office (project FIS2006-121 17). 
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